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ABSTRACT 

In  this  report,  the  heat  transfer  and  isothermal  flow  friction 
characteristics  for  seven  compact  heat  exchanger  surfaces  are  presented 
and  compared.   Five  of  the  surfaces  were  of  a  similar  triangular  fin 
configuration,  differing  in  material  and  method  of  construction,  and 
two  of  the  surfaces,  having  been  tested  previously,  were  of  a  modified 
rectangular  passage  configuration.   The  modified  rectangular  passage 
surfaces  showed  better  overall  performance  than  the  triangular  fin  sur- 
faces, and  matrices  constructed  from  perforated  material  showed  better 
overall  performance  than  similar  surfaces  constructed  from  non-perforat- 
ed material.   For  similar  matrix   configurations,  small  variations  in 
manufacturing  procedure  had  little  effect  upon  heat  transfer  and  flow 
friction  performance. 

The  heat  transfer  data  was  obtained  by  means  of  the  single-blow 
transient  testing  technique.   A  cylic  method  of  transient  testing  was 
also  investigated,  and  found  to  be  both  reliable  and  useful  in  that  the 
additional  method  of  testing  extends  the  Reynolds  number  range  for  the 
experimental  apparatus. 
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NOMENCLATURE 

English  Letter  Symbols 

A    Matrix  total  heat  transfer  area 

A,    Base  area  (plane  surface  area) 

A    Matrix  minimum  free  flow  area 
c 

A,.   Matrix  total  frontal  area 
fr 

A    Open  area  of  matrix  material  due  to 
perforations 

A    Matrix  solid  cross  sectional  area 
available  for  thermal  conduction 


\ 


H 


Matrix  solid  cross  sectional  conduction 
area  considering  effect  of  matrix  material 
perforations 


sq 

ft 

sq 

ft 

sq 

ft 

sq 

ft 

sq 

ft 

sq 

ft 

sq 

ft 

c  Specific  heat  (gas)  at  constant  Btu/(lbm  deg  F) 

p  pressure 

c  Matrix  material  specific  heat  Btu/(lbm  deg  F) 

s 

Cf  Flow  stream  capacity  rate  (mc  )  Btu/(hr  deg  F) 

C  Matrix  heat  capacity  (W  c  )  Btu/(deg  F) 

s  s  s 


Hydraulic  diameter  of  internal  ft 

passage 

E    Friction  power  expended  per  unit  of        hp/sq  ft 
surface  heat  transfer  area 

F    Frequency  of  sinusoidal  temperature        cycles/sec 
oscillation 

G    Matrix  flow  stream  mass  velocity  (m/Ac)    lbm/(hr  sq  ft) 


g    Proportionality  factor  in  Newton's         32.2(lbm  ft)^Lbf  sec  ) 
second  law 

h    Thermal  convection  surface  heat  transfer   Btu/(hr  sq  ft  deg  F) 
coefficient;  heat  transfer  power  per 
unit  area  per  degree  temperature  difference 

K    Loss  coefficient  for  sudden  flow  contrac-  dimensionless 
c   tion  at  matrix  entrance 


K  Loss  coefficient  for  sudden  flow 
g 

expansion  at  matrix  exit 

k  Unit  thermal  conductivity 

k  Matrix  material  thermal  conductivity 
s 

L  Total  matrix  flow  length 

m  Mass  flow  rate 

P  Pressure 

p  Matrix  porosity  (A  /A.  ) 

q  Heat  transfer  rate 

R  Universal  gas  constant  (53.35  for  air) 

r  Electrical  resistance 

r  Hydraulic  radius  (A  L/A.  ) 

t  Temperature 

u  Flow  velocity 

V  Matrix  volume 
m 

V  Matrix  material  volume  (corrected 
for  effects  of  perforations) 

W  Matrix  mass 
s 

x  Distance  from  matrix  inlet  in 
flow  direction 


dimensionless 

Btu/(hr  sq  ft  deg  F/ft) 

Btu/(hr  sq  ft  deg  F/ft) 

ft 

lbm/hr 

lbf/sq  ft 

dimensionless 

Btu/hr 

(ft  lbf)/(lbm  deg  R) 

ohm 

ft 

deg  F 

ft/sec 

cu  ft 

cu  ft 

lbm 
ft 


Greek  Letter  Symbols 


f*  Area  density;  compactness  (ratio  of  matrix  total 

heat  transfer  area  to  core  volume-A/V  )  sq  ft/cu  ft 

A  Difference  or  change  (time,  distance,  temperature) 

f[f  Fin  efficiency 

&  Time 

Ji  Fluid  viscosity  (dynamic) 

*/  Open  area  ratio  due  to  matrix  material 

perforations 

(P     Density 

*»      Ratio  of  free  flow  to  frontal  area  (A  /A,.  ) 

c  fr 

<P  Denotes  "function  of" 

Cdo  Angular  frequency  of  temperature  oscillation  rad/sec 


Sec,  hr. 
lbm/hr  ft 

lbm/cu  ft 
dimensionless 


Subscripts 

atm  Local  atmosphere 

ave  Average 

f  Fluid  (gas,  air) 

i  Initial,  inlet,  individual 

k  EqttitWent 

m  Mean,  matrix 

0  At  orifice 

s  Solid  (matrix  material) ,  static 

STD  Standard  (temperature  and  pressure) 

w  Wall  (solid  surface) 

x  Local  conditions 

1  At  inlet  (upstream  of  heating  elements) 

2  At  matrix  inlet  (downstream  of  heating) 

3  At  matrix  exit 
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Dimensionle8S  Groupings 

f        Fanning  friction  factor  (Ratio  of  wall  shear  stress  to  the 
fluid  dynamic  head) .   This  factor  is  plotted  versus  Nr  to 
illustrate  the  flow  friction  characteristics  of  a  matrix 
surface. 

ij 

j         Colburn  Factor  (  N$t  Npv-    )t  a  generalized  heat  transfer 

grouping.   This  factor  is  plotted  versus  Nr  to  illustrate  heat 
transfer  characteristics  of  a  matrix  surface 


Pr 
NR 
NSt 


Flow  parameter  (     ?/^C%co^ 


M 

N        Nusselt  number  (     /\c     )»  a  heat  transfer  modulus 

N 


Prandtl  number  (  *       /K  )  >  a  fluid  properties  modulus 
Reynolds  number  (  ^  /yU     ) ,  a  flow  modulus 
Stanton  number  (  /qcd    )»  a  heat  transfer  modulus 


h  A, 
Ntu       Number  of  heat  transfer  units  (   /#>•£-  )>  a  heat  transfer 

parameter 


R        Ratio  of  matrix  outlet  temperature  amplitude  to  matrix  inlet 
temperature  amplitude  during  cyclic  temperature  variation 

X       Longitudinal  heat  conduction  parameter  for  solid  matrix 
material  (  ****/*  Lcp      )• 

X«      Longitudinal  heat  conduction  parameter  for  perforated  matrix 
material  (  '<*  **/*/.*<*  >' 

T       Distance  parameter  (x/L) 

la  / 

«-      Time  parameter  (    /v/«  r       )  usea*  *n  Locke's  analysis 

<-      Time  parameter  (   y7_  )  used  in  Cyclic  technique 
&■>      Frequency  parameter  (   °  S£»  )  used  in  cyclic  technique 
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1.    INTRODUCTION 

A  compact  heat  exchanger  may  be  defined  as  a  heat  exchanger  having 
a  large  heat  transfer  surface  area  per  cubic  foot  of  core  volume. 
Arbitrarily,  heat  exchangers  with  more  than  200  square  feet  of  heat 
transfer  area  per  cubic  foot  of  core  volume  are  considered  compact. 
Applications  of  particular  interest  for  compact  heat  exchangers,  are  gas 
turbine  power  plants  in  aircraft,  automobiles,  and  small  seagoing  craft. 
In  all  these  applications  the  designer  is  limited  by  weight,  flow  frontal 
area,  core  volume,  and  pressure  drop  across  the  heat  exchange  r|<  An  effec- 
tive basis  for  presentation  and  comparison  of  heat  transfer  and  flow 
friction  data  for  various  heat  transfer  surfaces  and  core  configurations 
is  by  plotting  the  dimensionless  groupings-  Colburn's  heat  transfer  fac- 
tor, "j",  and  the  Fanning  friction  factor,  "f",  versus  Reynolds  number  for 
heat  exchanger  surfaces. 

The  data  presented  in  this  report  was  obtained  by  means  of  the 
single-blow  transient  testing  technique  described  in  the  two  sections 
following,  and  in  reference  [24],   The  objectives  of  this  thesis  were 
to  obtain  and  compare  flow  friction  and  heat  transfer  data  for  compact 
heat  exchanger  matrices  of  different  configuration,  to  correlate  this 
data  obtained  with  data  obtained  for  identical  surfaces  by  means  of 
a  steady  state  steam- to-air  testing  technique,  and  to  develop  a  testing 
technique  that  will  extend  the  range  of  transient  testing  to  higher 
Reynolds  numbers. 
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2.    SUMMARY  OF  THEORY 


A.   BACKGROUND 


Analytical  development  of  the  theory  for  the  single-blow  or 
transient  heating/cooling  problem  was  performed  in  the  1920' s  by  Anzelius 
(1926),  Nusselt  (1927),  Hausen  (1927-1929),  and  Schumann  (1929). 

There  was  some  duplication  of  effort  in  that  Nusselt  did  not  know 
of  Anzelius'  publication,  and  Schumann  was  unaware  of  the  work  of 
Anzelius  and  Nusselt.  Anzelius  and  Schumann  used  similar  procedures. 
Anzelius  represented  the  gas  and  solid  temperatures  as  unsolved,  re- 
latively simple  integrals,  where  Schumann  derived  and  solved  explicitly 
two  infinite  series  utilizing  Bessel  functions.   Schumann's  solation  was 
for  a  liquid  flowing  through  a  porous  solid,  however,  dimensional  analy- 
sis may  be  used  to  make  his  analysis  applicable  to  gasses  [24].  Hausen' s 
derivation  was  similar  to  Schumann's,  and  in  addition  Hausen  also  solved 
the  equations  for  cyclic  heating  and  cooling.   Nusselt 's  derivation  was 
more  simplified  than  Hausen *s  or  Schumann's,  but  it  is  worth  noting  that 
Nusselt  was  the  only  one  who  considered  the  problem  as  part  of  a  regenera- 
tor theory  [10] .  A  summary  of  the  single  blow  theory  as  it  applies  to 
a  gas  flowing  through  a  porous  solid  is  presented  below. 

B.   THEORY 

A  homogeneous,  porous  solid  has  a  fluid  flowing  through 
it.   Both  fluid  and  solid  are  at  the  same  constant,  uniform  temperature. 
After  a  step  change  in  the  entering  fluid  temperature,  the  temperature  of 
the  solid  and  the  fluid  are  to  be  determined  as  functions  of  time  and 
position  (along  the  flow  length  of  the  solid). 
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Starting  with  an  energy  balance  in  the  matrix: 


rncptp 


KSAS?£ 


WCp^'Hty) 


hb(tF-t3)dx 


-**fe*&*? 


Assumptions: 

(1)  The  fluid  properties  are  temperature  independent. 

(2)  The  fluid  flow  is  steady. 

(3)  The  porous  solid  is  homogeneous. 

(4)  The  thermal  conductivities  of  the  solid  and  the  fluid  are 
infinite  in  the  direction  perpendicular  to  the  fluid  flow. 

(5)  The  thermal  conductivity  of  the  fluid  is  zero  in  the  flow 
direction. 

Boundary  conditions: 

(1)  The  matrix  is  initially  at  uniform  temperature. 

(2)  At  time  zero,  there  is  a  step  change  in  the  entering  fluid 
temperature. 

(3)  The  matrix  boundaries  are  adiabatic. 
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Heat  rates  from  the  above  energy  diagram  are: 

(1)  Heat  absorbed  by  the  solid:  PsAa^S  ^«  dx 

(2)  Heat  transferred  to  the  solid  by  convection:  h  b  (tr  -tjjdx 

(3)  Heat  transferred  from  the  fluid:  |Y|  Cp  ^^d 


Sx 


K 


(4)  Heat  transferred  in  the  solid  by  conduction:   -  K5A5  "-^^Jx 
The  net  resulting  energy  balances  for  fluid  and  solid  are: 

(1)  Fluid:   mCplj-^cU  +hfc  (tp-ts)dx  -  O  (2-1) 

(2)  Solid:      PsAsCfcH^x  -  k3As^^hL(tr-tis)j<(2.2) 
Now  introduce  the   following  dimensionless  parameters: 

(1)  Time  parameter 

T-  %k  (•  -  *  t) 

Where: 

r\  =  unit  conductance  for  convective  heat  transfer 


A  =  matrix  heat  transfer  area  (sq  ft) 
V$C^  ■  matrix  capacity  (Btu/deg  F) 
C/  ■  time  (hours) 
V/f  =  mass  of  fluid  contained  within  matrix  (lb  m) 
m  =  fluid  mass  flow  rate  (lbm/hr) 

*   ■  distance  from  matrix  inlet  in  flow  direction  (ft) 
L  =  matrix  flow  length  (ft) 
Introducing  the  fluid  specific  heat,  and  regrouping,  (2-3)  becomes: 

T=M-s-M-f-  wt^  (2.4) 

V/rCp 
The  term    L   ■    is  very  small  since  the  thermal  capacity  of  the  matrix 

is  much  greater  than  the  thermal  capacity  of  the  fluid  contained  within 

it.   Thus,  the  dimensionless  time  parameter  can  be  reduced  to: 
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T  s  v£T5  e  <2-5> 

(2)  Position  parameter; 

Since  the  Number  of  Transfer  Units,  NtU  s  fi*ci»  : 

Z*  Kfe.t\-f-  (2-7) 

(3)  Conduction  parameter: 

v    ksAs 

A  *   .  ■   ;  (2-8) 

Where: 

k  =  matrix  thermal  conductivity  (Btu/hr  sq  ft  deg  F/ft) 

s 

A  ■  solid  matrix  cross  sectional  area  available  for  thermal 
s 

conduction  (sq  ft) 

Introducing  the  above  dimensionless  parameters  into  the  heat  balance 

equations  for  the  fluid  and  solid,  and  rearranging  yields: 

Fluid:   ^=ts-tf  (2-9) 

Solid:   SjfcA.  tf  .is+XNi£t<-|^  (2-10) 

If  thermal  conduction  in  the  solid  is  assumed  zero  in  the  direction 

of  flow,  the  above  equations  simplify: 

|*f-  U~tp  (2-1D 

¥?-u-u  (2-12) 

From  equations  (2-11)  and  (2-12),  Schumann,  utilizing  Bessel  func- 
tions, developed  his  temperature-time-distance  relationships: 

tf-t;  &      Tt^f (2"13) 

tF.-t;  '-'-  &  >  Z    dC*l]" 

Plots  of  these  functions  versus  1   appear  in  Figure  1.  Using  these 


][ 


curves  to  determine  Ntu  from  experimental  data  would  be  a  tedious  task 
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involving  infinite  series  computations  and  curve  comparisons.  However, 
in  his  maximum  slope  technique,  Locke  [15]  differentiated  the  theoretical 
curves,  plotting  slope  versus  ^/Ntu  for  various  values  of  Ntu.   Locke's 
expression  for  slope  is: 


Ntu 


±n±\ 


(2-15) 


The  downstream  fluid  temperature  is  taken  at  x  =  L  so  that  21  be- 
comes Ntu  and  tf  =  tf  .   Locke  also  evaluated  the  maximum  slope  of  the 
temperature  curve  for  each  Ntu,  and  plotted  maximum  slope  versus  Ntu,  a 
single  curve.  Figures  2  and  3  show  Locke's  curves. 

Both  Schumann  and  Locke  assumed  thermal  conductivity  of  the  solid, 
in  the  flow  direction,  to  be  zero.  The  effect  of  longitudinal  thermal 
conductivity  is  small  in  the  high  Reynolds  number  range,  but  it  is  con- 
siderable and  should  be  taken  into  account  in  the  low  Reynolds  range. 
To  account  for  longitudinal  thermal  conduction,  equations  (2-9)  and 
(2-10)  must  replace  equations  (2-11)  and  (2-12).  Thus,  the  maximum 
slope  will  depend  upon  both  Ntu  and  X  .  Howard  [9],  utilizing  a 
digital  computer  and  finite  difference  technique,  tabulated  and  plotted 
Ntu  versus  maximum  slope  for  various  values  of  A  .  Howard's  curves 
appear  in  Figures  4  through  6.   By  means  of  a  properly  designed  experi- 
mental apparatus  which  satisfies  Howard's  idealizations  and  boundary 
conditions,  the  maximum  slope  of  the  heating  or  cooling  curve  and  the 
conduction  parameter  can  be  determined.   Ntu  can  then  be  found  using 
Howard '8  data  shown  in  Figure  5. 

In  the  vicinity  of  Ntu  ■  2,  it  can  be  seen  from  Figure  6  that  there 
can  be  considerable  error  in  Ntu  for  a  small  error  in  maximum  slope. 
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In  addition,  it  is  more  difficult  to  accurately  determine  maximum  slope 
in  the  high  Reynolds  number  region,  as  the  maximum  slope  occurs  earlier 
in  time  as  Ntu  decreases,  and  the  determination  of  the  maximum  slope 
becomes  difficult,  having  large  uncertainties  due  to  the  inherent  time 
response  of  the  instrumentation. 

Another  theory  for  determining  the  heat  transfer  characteristics  of 
a  porous  solid,  based  upon  a  sinusoidal ly  varying  fluid  temperature  at 
the  matrix  inlet,  was  developed  by  Bell  and  Katz  [2].   Starting  as  be- 
fore with  an  energy  balance  in  the  matrix: 


mCpt-F 


tefa+^^W*) 


hb(fcF-ts>k 


Assumptions: 

(1)  All  cross  sections  of  the  solid  normal  to  the  flow  direction 
are  uniform. 

(2)  The  fluid  velocity  is  uniform  throughout  the  matrix. 

(3)  The  fluid  properties  are  temperature  independent. 

(4)  Thermal  conductivity  of  the  solid  and  fluid  is  infinite  in  the 
direction  normal  to  fluid  flow. 
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(5)  Thermal  conductivity  of  the  fluid  and  solid  is  zero  in  the 
direction  parallel  to  the  fluid  flow. 

(6)  The  matrix  boundaries  are  adiabatic. 

(7)  The  convective  heat  transfer  coefficient  remains  constant. 
Heat  rates  from  the  above  heat  balance  are: 

(1)  Heat  absorbed  by  the  solid:  Ps /4$  C5  ^=^ 

(2)  Heat  transferred  to  the  solid  by  convection:    hb  ^£f  -  £.3)  <dx 

(3)  Heat  absorbed  by  the  fluid:  fftf pfejs£/v   2 £ p  ,  \ 

(4)  Heat  transferred  in  the  solid  by  conduction:         q 
The  net  resulting  energy  balances  for  fluid  and  solid  are: 

(1)  Fluid:  rinCp  (l|£jlt+!^e)«/)b(fcF-£9)<Jx     (2-16) 

(2)  Solid:      psA5CS^K    =hb(fcf-t5)dx  <2-17) 
Inlet  and   initial  conditions  are: 

(1)  Inlet:       t  f  (Oj&)    «   £  m  +  A  **t\    OJ*&  (2-18) 

(2)  Initial:  £s  (p)         *   t5o  (2-19) 

(3)  Initial:  £p(xjO)-  £o(*)  (2"20) 
where  A  is  the  amplitude  of  the  sinusoidally  varying  temperature. 

Now  introduce  the  following  dimensionless  parameters: 

(1)  Time:   c     ^  -  ^j*  (2-21) 

(2)  Distance:     t  "  *L  (2-22) 

(3)  Frequency:    60  *.  °>°L  (2-23) 

(4)  Flow:        fi  *   "^£f  (2-24) 

WsCsTOe 

(5)  Heat  transfer:  ...      U  (2-25) 

Where: 

U  =  flow  velocity  (ft/sec) 

©  =  time  (sec) 

L  =  matrix  flow  length  (ft) 
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X  ■  distance  from  matrix  inlet  in  flow  direction  (ft) 

60o=  angular  frequency  of  temperature  oscillation  (rad/sec) 
Introducing  the  above  dimensionless  parameters  into  equations  (2-16)  through 
(2-20)  and  rearranging  yields: 

(1)  §f*+  &*  +  hltli(t.F-t3)-0  (2-26) 

(2)  W-M<oMt«<:*r-i5)  =  o  (2-27) 

(3)  LrtO,*)  ~   tto+LHh   Cu^  (2-28) 

(*)  ts(O)    =  t5.  (2-29) 

(5)  tf  (f,0)  -  ±o(f)  (2-30) 
This  system  has  the  periodic  solution: 

tf(^!)s^4Ae    s/h  fdo*E-jffT)  (2-31> 

t s  ( f  >  X)  =  t  hn  +  A°t  £    [n  Kl£u  S.'n  <*j«  -jSfJ  -  COS  (So^  -^f)7   (2-32) 
Where:  ^ 


A  simple  quantity  to  determine  experimentally  is  RA,  the  ratio 


~A] 
of  the  amplitude  of  outlet  and  inlet  gas  temperatures.   Straightforward 

calculation  shows: 


-ct 


Nln 


Ka  "    &  (2-33) 

A  plot  of  Ntu  versus  R.  for  various  values  of  M  is  shown  in  Figure  7. 
The  same  experimental  apparatus  can  be  used  for  this  cyclic  operation 
as  is  used  for  the  maximum  slope  method.  The  only  addition  needed  is 
a  method  to  sinusoidally  vary  the  temperature  of  the  air  entering  the 
matrix. 

Thus,  heat  transfer  data  can  be  determined  for  a  wide  range  of 
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Reynolds  numbers  utilizing  two  testing  methods  on  a  single  experimental 
apparatus;  the  single-blow  technique  producing  heat  transfer  data  for 
low  Reynolds  numbers,  and  the  cyclic  technique  producing  data  for  high 
Reynolds  numbers. 


22 


3.    EXPERIMENTAL  TECHNIQUES 

A.   Single  Blow  Technique 

For  Reynold's  numbers  up  to  approximately  1,000,  the  results 
of  Howard's  work  [9]  were  used  in  determining  the  heat  transfer  data, 
Ntu  and  j,  appearing  in  this  report.   In  order  to  utilize  Howard's 
theoretical  data,  the  experimental  apparatus  had  to  conform  to  the  ideali- 
zations stated  in  reference  [9]. 

Idealizations: 

(1)  The  fluid  flowing  through  the  matrix  remains  steady  and 
uniform  in  velocity  and  temperature  at  any  cross  section.   The  matrix 
thermal  conductivity  is  infinite  in  the  direction  normal  to  fluid  flow, 
and  finite  in  the  direction  parallel  to  fluid  flow.   Thus  the  problem  is 
one -dimensional  in  space. 

(2)  The  thermal  capacity  of  the  matrix  is  large  compared  with 
the  fluid  contained  within  it.   This  restricts  the  fluid  to  a  gas,  and 
means  that  there  will  be  no  time  dependent  terms  in  the  equations  for 
the  fluid. 

(3)  The  thermal  properties  of  the  fluid  and  matrix  are  cons- 
tant and  uniform. 

(4)  The  convective  heat  transfer  coefficient  is  some  suitable 
average  and  remains  constant. 

(5)  At  time  zero,  the  change  in  fluid  temperature  is  a  step 
change,  with  the  matrix  and  its  entrained  fluid  initially  at  uniform 
and  equal  temperatures. 

The  test  apparatus  utilized  for  this  report  satisfies  the  above 
idealizations  as  described  below: 
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(1)  A  square  cross  section  contracting  cone  and  screen  type 
flow  straighteners  provided  the  necessary  uniform  velocity  profile.   In 
addition,  the  matrix  test  section  was  located  as  far  upstream  as  possible 
to  take  advantage  of  this  uniform  velocity  profile.  The  flow  rate 
measuring  apparatus  was  located  downstream.  The  heater  configuration, 
Figure  10,  insured  a  uniform  temperature  profile  due  to  even  distribu- 
tion of  heat  over  the  flow  cross  section.  The  thermal  conductivity  in 
the  direction  parallel  to  fluid  flow  is  taken  care  of  by  the  use  of  the 
conduction  parameter  X. 

(2)  The  use  of  a  gas  as  the  fluid  flowing  satisfies  idealiza- 
tion (2). 

(3,4)  A  small  temperature  variation  will  insure  uniform  solid 
and  fluid  physical  properties.  At  80°F,  a  temperature  variation  of  + 
10°F  will  result  in  an  air  viscosity  variation  of  approximately  +  1.5%. 
The  changes  in  specific  heats  and  thermal  conductivities  are  negligible 
for  a  +  10°F  change  in  temperature. 

(5)  Various  methods  have  been  devised  to  produce  a  step  change 
in  matrix  inlet  temperature.  The  electric  heater  arrangement  used  in 
this  report  appears  to  produce  the  closest  approximation  to  the  desired 
step  change  with  a  minimum  effect  on  flow  properties  [24]. 

A  schematic  diagram  of  the  experimental  apparatus  is  shown  in 
Figure  8.   Photographs  of  various  components  are  shown  in  Figures  9 
through  14.   The  system  is  comprised  of  an  entrance  section  with  a  contract- 
ing cone  and  a  screen  flow  straightener,  a  heat  section  consisting  of  28 
resistance  type  wire  heaters  made  of  nichrome  wire  in  grid  fashion,  a 
matrix  test  section,  an  orifice  metering  section,  a  flow  control  valve 
and  prime  mover,  and  pressure  and  temperature  measuring  and  recording 
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systems.   Each  of  these  components  are  further  described  in  Appendix 
A.   The  experimental  procedure  for  the  single  blow  technique  is  as 
follows: 

Air  is  drawn  through  the  system  by  means  of  a  30  horsepower  Spencer 
Turboblower.   The  heater  section  is  designed  to  heat  the  air  to  approxi- 
mately twenty  degrees  above  ambient  temperature.  The  mass  flow  rate  is 
measured  by  means  of  an  ASME  orifice  meter  with  removable  orifice  plates. 

Flow  friction  data  is  obtained  by  measuring  the  isothermal  pressure 
drop  across  the  matrix  and  the  mass  flow  rate  through  the  orifice  meter. 

Heat  transfer  data  is  obtained  as  follows:   the  heaters  are  energi- 
zed to  heat  the  air  to  approximately  twenty  degrees  above  ambient  tempera- 
ture. When  the  matrix  has  been  heated  to  a  uniform  temperature,  the 
heaters  are  quickly  de-energized  and  the  time-temperature  history  of  the 
air  leaving  the  matrix  is  recorded  on  a  Minneapolis  Honeywell  "Brown" 
strip  chart  recorder. 

The  following  data  is  recorded  for  each  test  run: 

(1)  Atmospheric  pressure  Patm 

(2)  Static  pressure  upstream  of  matrix  P 

(3)  Pressure  drop  across  the  matrix  A  P 

m 

(4)  Static  pressure  upstream  of  orifice  P 

(5)  Pressure  drop  across  the  orifice  A  P 

(6)  Temperature  upstream  of  orifice  t 

(7)  Orifice  diameter  d 

o 

(8)  Time- temperature  cooling  curve 

(9)  Recorder  chart  speed 
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B.   Cyclic  Technique: 

To  obtain  heat  transfer  data  by  means  of  the  cyclic  tempera- 
ture variation  method,  the  same  testing  apparatus  can  be  used,  the  only 
modification  being  in  the  power  input  to  the  heaters.   The  idealizations 
set  forth  by  reference  [2]  are  the  same  as  those  imposed  by  reference  [9] 
and  shown  on  the  previous  pages  with  two  exceptions: 

(a)  The  thermal  conductivity  of  the  solid  is  zero  in  the 
direction  parallel  to  the  fluid  flow. 

(b)  The  temperature  of  the  matrix  inlet  varies  sinusoidal - 
ly.   To  satisfy  idealization  (a),  this  testing  method  is  restricted  to 

the  high  Reynolds  number  ranges  (  [\/^  >  500),  where  Ntu  is  small  (Ntu  ^,  5), 
and  the  effect  of  longitudinal  thermal  conduction  is  small  and  can  be 
neglected.   (See  Figure  15  and  Chapter  5). 

To  satisfy  idealization  (b) ,  sinusoidally  varying  power  is  ap- 
plied to  the  heaters  by  means  of  a  1.5  KW  Amplidyne  receiving  field  excita- 
tion from  a  variable  frequency  function  generator.   This  system  is  des- 
cribed further  in  Appendix  B. 

To  obtain  heat  transfer  data,  the  temperatures  of  the  matrix  inlet 
and  exit  are  recorded  on  the  strip  chart  recorder.   The  amplitudes  of 
the  cyclic  variation  of  these  two  temperature  curves  are  then  measured 
and  Ntu  is  determined  from  the  ratio  of  the  amplitudes  and  the  dimension- 
less  parameter,  M.   Figure  15  shows  a  sample  chart  trace  for  a  cyclic 
test. 

The  same  data  are  recorded  as  for  the  maximum  slope  technique. 

Data  reduction  relationships  for  both  of  these  transient  testing  tech- 
niques are  described  in  Appendix  B. 
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4.    DESCRIPTION  OF  MATRICES 

The  surfaces  investigated  in  this  report  are  of  two  major  varieties 
modified  rectangular  passage  and  triangular  fin. 

The  modified  rectangular  passage  surfaces,  previously  examined  by 
Bannon  [1],  and  Piersall  [22],  are  constructed  of  formed  sheets  separat- 
ed by  flat  splitter  plates  as  shown  below.   The  matrix  approximates  one 
made  up  of  rectangular  passages.   The  rectangular  passages  have  an  aspect 
ratio,  OC  ,  of  approximately  7: 


One  parallel  plate  matrix  is  made  from  .0022  inch  thick  perforated 
nickel  and  the  other  from  .002  inch  thick  solid  nickel.   The  perforations 
are  of  type  160/40  TV  and  are  described  further  in  Appendix  C. 

The  five  triangular  fin  matrices,  provided  by  Solar,  A  Division  of 
International  Harvester,  San  Diego,  were  of  varied  material  and  method 
of  construction.   Geometric  and  physical  data  for  these  matrices  appear 
in  Figures  16  through  22.   Two  were  made  from  nickel;  one  had  fins  made 
from  .005  inch  thick  solid  material,  the  other  had  fins  made  from  type 
80/20  T  perforated  nickel  material  .0047  inches  thick.   The  other  three 
triangular  fin  matrices  were  made  from  .005  inch  thick  type  430  stain- 
less steel,  and  had  no  perforations.   Of  these  stainless  steel  matrices, 
one  had  fins  with  very  sharp  bends,  (all  others  had  approximately  sinusoi- 
dal corrugations) ,  one  had  maximum  fillet  braze  joints  where  fins  and 
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splitters  touched,  (all  others  were  not  bonded  by  any  means),  and  the 
last  stainless  steel  surface  was  similar  in  shape  and  construction  to 
its  solid  nickel  finned  counterpart.  All  of  the  triangular  fin  matrices 
had  a  fin  height  of  .1  inches  and  10  fins  per  inch.  All  of  the  matrices 
measured  3-1/16  inch  high  X  3-1/16  inch  wide.   The  modified  rectangular 
passage  matrices  were  2  inches  in  flow  length,  while  the  triangular 
finned  ones  were  3  inches  long.   Geometric  and  physical  data  for  all 
matrices  appear  in  Figures  16  through  22. 
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5.    PRESENTATION  OF  RESULTS 

Heat  transfer  and  flow  friction  data  for  all  surfaces  tested  are 
shown  in  TABLES  I  through  VII  and  Figures  23  through  35.   In  Figures 
23  through  31  the  Colburn  heat  transfer  modulus,  j,  and  Fanning  fric- 
tion factor,  f,  are  plotted  versus  Reynolds  number.   The  Reynolds  number 
was  based  on  hydraulic  diameter.   In  Figure  32,  all  surfaces  are  compared 
by  a  plot  of  j/f  versus  Reynolds  number,  and  Figure  33  is  a  plot,  for 
all  surfaces,  of  heat  transfer  power  versus  flow  friction  power  on  a 
unit  area  basis. 

Figure  34  is  a  plot  of  Colburn  j  versus  Reynolds  number  for  the  brazed 
stainless  steel  triangular  fin  surface.  The  j  values  for  this  plot  were 
obtained  by  means  of  both  the  single-blow  technique,  and  the  cyclic  tempera- 
ture variation  technique  described  in  Chapters  2  and  3.   For  N  <  1000 
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the  single-blow  technique  was  employed,  and  for  N  >  500  the  cyclic  tech- 
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nique  was  used.   This  was  the  only  surface  tested  by  means  of  the  cycli- 
cal temperature  variation  technique, 

Figure  15  shows  a  plot  of  Colburn  j  versus  Reynolds  number  for  the 
perforated  nickel  rectangular  passage  matrix.   This  plot  shows  the  ef- 
fects of  varying  the  longitudinal  conduction  parameter,  X  .   One  curve 
results  from  assuming  A  =  0,  the  second  curve  results  from  the 
calculated  in  reference  [22]  and  described  in  Appendix  C,  and  the  third 
curve  results  from  an  electrical  analogy  method  calculation  of  A  which 
is  also  described  in  Appendix  C. 

Figure  35  is  a  plot  of  Colburn  j  versus  Reynolds  number  for  the 
perforated  nickel  triangular  fin  matrice.   the  data  for  this  curve  was 
obtained  by  means  of  both  the  single  blow  technique  and  a  steady  state 
steam  to  air  technique  described  in  reference  [29]. 
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In  Figure  36,  experimental  and  theoretical  data  are  compared  for 
the  brazed  triangular  fin  matrix. 
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6.    EXPERIMENTAL  UNCERTAINTIES 

The  experimental  uncertainties  for  the  heat  transfer  and  friction 
data  appearing  for  this  report  were  determined  by  the  method  described  in 
reference  [13]  and  shown  in  Appendix  D.   These  uncertainties  arise  from 
possible  errors  in: 

(1)  Physical  constants  -  These  physical  constants  were 
obtained  from  references  [6],  [7],  and  [8],  and  their  uncertainties  are 
estimated  to  be  as  follows: 

ks      ±  -5% 

C5    ±  .5% 

CP   +  .5% 

Npr  ±  2*07« 
yu    +  1.0% 

(2)  Geometric  measurements  -  These  errors  are  due  to  con- 
struction inaccuracies  and  errors  in  linear  measurement  and  are  esti- 
mated to  be  as  follows: 

Ac,Afr,A,A&     ±l-« 

L   ±  .5% 

W$   negligible 

(3)  Instrumentation  -  Assuming  adequate  manufacturer's  cali- 
bration a  manometers  and  thermocouple  wire,  instrumentation  errors  will 
be  assumed  to  be  due  to  instrument  fluctuations  (manometers) ,  and  size 
of  scale  divisions  (manometers  and  potentiometer). 

Pressures: 

Ps    +  1.07. 
APm        ±   1-07. 
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p0    +  1.07. 

AP0   ±  l-w 

Path        negligible 
The  error  for  the  manometers  depends  upon  the  range  of  the  pressure 
being  read.   In  the  range  1  to  3  inches  of  water,  the  error  is  .17%  to 
.57.  in  the  range  3  to  5  inches  of  water,  the  error  is  1.0%  to  1.7%. 
Temperatures: 

to  ±  L5  «/«g  F 
Maximum  Slope; 

The  measurement  of  maximum  slope  involves  the  drawing  of 
a  tangent  and  the  linear  measurement  of  three  quantities  X  ,  y  ,  and 
inn   ,  as  described  in  Appendix  A.  The  percent  error  in  maximum  slope 
is  estimated  as  +  2.0%.   Entering  Figure  (7)  with  this  value  at  Ntu  =  3.0 
and  >v  =  0,  obtain  Ntu  error  as  +  7.2%.  At  Ntu  -  25  and  X  -  .05,  Ntu 
error  is  +  10.0%  for  a  2.07»  error  in  maximum  slope. 

A  summary  of  extreme  values  of  uncertainty  intervals  calculated  by 
the  method  shown  in  Appendix  D  is  given  below: 

*Y|    +  1.0% 
|s/r        ±  2-37- 

j       ±  I0* 
f      ±  47» 


32 


7.    DISCUSSION  OF  RESULTS 

Heat  transfer  and  flow  friction  results  appear  in  Figures  23  through 
36,  with  comparisons  of  all  surfaces  in  Figures  31  through  33. 

In  general,  the  triangular  fin  surfaces  all  had  similar  j  and  f 
curves;  the  nickel  matrices  having  higher  j  and  j/f  values  than  the  stain- 
less steel  matrices.  The  high  j  and  j/f  for  the  nickel  surfaces  are  due 
chiefly  to  the  high  thermal  conductivity  of  nickel  -  36.0  compared  to 
12.8  for  430  stainless  steel  (at  80°F).   The  j,f,  and  j/f  characteristics 
for  the  triangular  fin  surfaces  are  in  agreement  with  the  theory  presented 
in  reference  f 16 ]  for  laminar  flow  through  idealized  cross  sections  (see 
Figure  36). 

The  three  stainless  steel  cores  were  quite  similar  in  cross  section, 
however,  one  might  expect  to  find  small  differences  in  their  heat  trans- 
fer and  flow  friction  characteristics  due  to  slight  differences  in  geo- 
metry. The  rounded  corners  produced  by  the  sine  wave  corrugations  in 
two  of  the  surfaces  should  reduce  the  effective  surface  area  for  heat 
transfer  and  wall  friction  as  mentioned  in  reference  [16].  The  experi- 
mental results  seem  to  agree  in  part  with  the  above  -  the  sharp  cornered 
triangular  passage  geometry  yielded  higher  friction  factors  than  the 
sinusoidal  fin  geometry,  with  the  brazed  sinusoidal  fin  core  producing 
intermediate  friction  factors.   The  heat  transfer  characteristics,  how- 
ever, showed  the  j  curve  for  the  brazed  core  slightly  higher  than  the  j 
curves  for  the  unbrazed  ones.  There  are  two  reasons  for  this:  First  of 
all,  the  heat  transfer  area  calculations  for  the  brazed  core  took  into 
account  the  loss  of  metal  surface  area  due  to  brazing.   In  the  unbrazed 
sinusoidal  fin  core,  a  portion  of  surface  area  similar  to  the  "lost  area" 
was  relatively  ineffective  for  heat  transfer,  but  not  accounted  for  in 
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calculating  j.   Secondly,  the  fin  shape  and  matrix  construction  for  the 
brazed  core  actually  produced  a  geometry  that  is  a  cross  between  a  tri- 
angular passage  and  a  rectangular  passage,  with  a  rectangular  passage 
having  higher  values  of  j  and  j/f  than  a  triangular  passage  as  shown  in 
reference  [16]. 

It  should  be  noted  that  these  three  stainless  steel  matrices  yielded 
heat  transfer  and  flow  friction  results  that  differed  only  slightly,  in- 
dicating that  small  deviations  in  manufacturing  procedure  may  have  little 
effect  upon  heat  transfer  and  flow  friction  performance. 

The  rectangular  passage  cores  had  heat  transfer  and  friction  charact- 
eristics similar  to  the  triangular  fin  cores;  however,  their  heat  trans- 
fer curves  (j  versus  N  )  were  less  steep,  the  perforated  material  pro- 
ducing a  markedly  higher  curve  than  the  solid  material.   The  j/f  flow 

area  "goodness"  factors  for  these  surfaces  increased  as  N   increased, 

R 

while  j/f  for  the  triangular  fin  surfaces  decreased  with  increasing 
Reynolds  number  (See  Figure  32).   On  the  h    versus  E    plot,  Figure 
33,  rectangular  passage  surfaces  showed  definite  superiority  over  the 
triangular  fin  surfaces,  with  the  perforated  material  producing  the  best 
results. 

Figure  15  demonstrates  the  influence  of  conduction  parameter  upon 
the  heat  transfer  modulus,  j.   Curves  such  as  those  shown  dotted  can 
result  from  an  improper  assumption  (  X  =  0  in  low  Reynolds  range) ,  ap- 
proximations for  X  (Appendix  C) ,  or  an  inaccurate  value  for  k  . 

s 

Referring  to  Figure  7,  it  can  be  seen  that  there  are  areas  where  a 
small  error  in  Amplitude  Ratio,  R.,  would  result  in  a  large  correspond- 
ing error  in  Ntu.   However,  the  frequency  of  temperature  oscillation, 
60©  »  can  be  chosen  so  as  to  produce  a  flow  parameter,  M,  (M  =  mc  /£^Cs&^) 
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that  will  result  in  Ntu  being  determined  with  minimum  error  for 
a  given  error  in  R  .   For  instance,  at  Ntu  ■  1.5,  C)Q should  be  chosen  so 
that  M  is  less  than  .35,  since  values  of  M  near  .7  in  this  case  result 
in  a  high  uncertainty  of  Ntu.  This  flexibility  is  not  inherent  in  the 
single-blow  technique  (where  Ntu  «  2) . 

Figure  35  shows  heat  transfer  data  for  the  brazed  triangular  fin 
matrix.   Both  the  single-blow  and  cyclic  techniques  were  employed,  thus 
producing  data  up  to  a  Reynolds  number  of  1850  with  an  area  of  overlap 
between  Reynoldd  numbers  380  and  560.   There  seems  to  be  fairly  good 
correlation  between  the  two  testing  methods  with  the  exception  of  the 
data  points  at  Reynolds  numbers  450  and  560.   These  points  were  the  re- 
sult of  employing  the  single-blow  technique  in  the  area  of  high  N  and 
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low  Ntu  (1.65  and  1.44),  where  Ntu  is  subject  to  error  as  explained  in 
Chapter  2.   The  data  from  the  cyclic  technique  was  also  subject  to  some 
error  due  chiefly  to  the  cyclic  power  input  to  the  heaters.  As  can  be 
seen  on  the  sample  chart  trace  in  Figure  37,  there  was  some  distortion 
in  the  temperature  waveforms  making  them  slightly  different  than  true 
sinusoids.   This  distortion  is  due  to  inaccurate  signal  reproduction  in 
the  amplidyne.   There  was  also  a  problem  in  instrumentation  -  the  tempera- 
ture waveforms  were  recorded  separately  by  means  of  the  same  instrument 
instead  of  simultaneously.   Bell  and  Katz  [2]  utilized  a  dual  channel 
instrument  recording  both  temperatures  simultaneously,  reversing  the  in- 
put connections  midway  through  each  run  so  as  to  record  each  temperature 
on  both  channels,  thus  eliminating  the  effects  of  unequal  channel  sensi- 
tivities.  The  thermocouple  configuration  described  in  Appendix  B  does 
not  allow  for  simultaneous  temperature  measurement  and  recording. 
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In  Figure  35  data  from  the  single-blow  technique  is  plotted  with 
data  from  a  steady  state  stream  to  air  technique.   In  both  tests  the 
identical  type  80/20T  perforated  nickel  triangular  fin  surface  was  in- 
vestigated.  The  steam  to  air  test  was  conducted  on  a  cross  flow  heat 
exchanger  with  the  surface  to  be  tested  on  the  air  side.   The  core  con- 
figuration for  the  single-blow  test  is  described  in  Chapter  4  and  Figure 
17.   Theoclitus  in  reference  [28]  utilized  a  dual  testing  technique  similar 
to  this  in  determining  fin  efficiency,  7/f  ,  and  total  surface  temperature 
effectiveness  /(0  ,  for  various  surfaces.   Theoclitus  found  for  simple 
finned  surfaces,  that  the  theories  for  fin  efficiency  and  total  surface 
temperature  effectiveness  are  reliable  since  he  obtained  a  close  corre- 
lation of  data  from  the  two  methods  of  testing.   For  the  steady  state 
technique,  theoretical  values  of  /f^-and   ff^   are  used  in  determining  the 
heat  transfer  characteristics,  Ntu  and  j,  while  [l~   and  ifo  have  values 
of  unity  when  the  transient  technique  is  employed,  due  to  lack  of  a 
temperature  gradient  in  the  fin  in  the  direction  perpendicular  to  fluid 
flow  (See  Chapter  3). 

Figure  35  indicates  a  fairly  good  correlation  of  data  for  the  two 
testing  methods. 

Conclusions: 

(1)  For  matrices  of  the  same  geometry  and  material,  small  devia- 
tions in  manufacturing  procedure  have  little  effect  upon  overall  heat 
transfer  and  friction  performance. 

(2)  The  single-blow  and  cyclic  transient  testing  techniques 
utilized  for  this  report  are  reliable  in  that  the  experimental  results 
obtained  agree  with  results  obtained  by  means  of  steady  state  testing 
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and  with  theory. 

(3)  The  single-blow  and  cyclic  techniques  supplement  each  other 
in  producing  heat  transfer  data  over  a  wide  range  of  Reynolds  numbers. 
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8.    RECOMMENDATIONS  FOR  FURTHER  WORK 

Refinement  of  the  cyclic  method  of  testing  is  needed  both  in  theory 
and  in  testing  apparatus.  The  theory  could  be  extended  to  include  the 
effect  of  longitudinal  thermal  conduction,  this  involving  the  use  of 
another  non-dimensional  parameter  in  place  of  M,  and  a  theoretical  solu- 
tion using  methods  similar  to  Howard's  [9]  or  Moreland's  [20]. 

There  are  a  number  of  areas  for  improvement  in  the  test  apparatus: 

(1)  A  modification  is  needed  if  higher  flow  rates  are  to  be 
employed.   Either  a  larger  pipe  diameter  at  the  orifice 
meter,  or  another  means  of  measuring  flow  rate  is  needed 
to  increase  the  flow  rates  employed. 

(2)  A  more  reliable  power  source  is  needed  for  the  heaters  so 
that  a  truer  sinusoidal  temperature  variation  results. 

(3)  The  temperature  measuring  system  could  be  modified  so 

that  the  matrix  inlet  and  exit  temperatures  can  be  recorded 

simultaneously. 
It  is  also  recommended  that  more  surfaces  be  made  available  for 
testing  by  both  transient  and  steady  state  methods  so  that  data  is  avail- 
able for  a  wider  range  of  Reynolds  numbers,  and  the  reliability  of  the 
two  methods  may  be  compared. 
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Figure  5.   Ntu  versus  Maximum  Slope 
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Figure  7.   Ntu  versus  Amplitude  Ratio 
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Figure  15.   Effects  of  Longitudinal  Conduction  Parameter 
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Matrix  Material 

Specific  Heat,  c      > 

Thermal  Conductivity,  k 

Material  Thickness 

Total  Heat  Transfer  Area 

Frontal  Area,  A_ 

f  r 
Conductive  Area,  A 

Free  Flow  Area,  A 

Matrix  Flow  Length, 

Matrix  Volume,  V 

Matrix  Density, 

Compactness, 

Porosity,  p 

Hydraulic  Diameter 


m 


m 


m 


SOLAR  #1 

Nickel 

(Btu/lbm  deg  F) 

.106 

(Btu/hr  ft  deg  F) 

36.0 

(in) 

.005 

,A  (sq  ft) 

.  12.768 

(sq  ft) 

.07022 

(sq  ft) 

.01078 

(sq  'ft) 

.05944 

(ft) 

.24917 

(cu  ft) 

.0175 

(lbm/cu  ft) 

86.1378 

(sq  ft/cu  ft) 

729.7026 

.8464 

(ft) 

.0046395 

Figure  16  Geometric  and  Physical  Properties 
for  Solar  #1 
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Matrix  Material 

Specific  Heat,  c 

Thermal  Conductivity,  k 

Material  Thickness 

Total  Heat  Transfer  Area 

Frontal  Area,  A, 

f  r 
Conductive  AreaB  A 

s 
Free  Flow  Area,   A 

Matrix  Flow  Length ,  L 

Matrix  Volume,  V 

Matxix  Density, 

Compactness 9 

Porosity,  p 

Hydraulic  Diameter 

Perforated  Fin  Material 


,A 


m 


m 


SOLAR  #2 

Nickel 

(Btu/lbm  deg  F) 

.106 

(Btu/hr  ft  deg  F) 

36.0 

(in) 

.005 

(sq  ft) 

12.093 

(sq  ft) 

.07022 

(sq  ft) 

.01034 

(sq  ft) 

.05987 

(ft) 

.25583 

(cu  ft) 

.01797 

(lbm/cu  ft) 

70.3099 

(sq  ft/cu  ft) 

673.143 

.85281 

(ft)' 

.0046272 

.0047"  thick 

80/20T 

Figure  17  Geometric  and  Physical  Properties 
for  Solar  #2 
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Matrix  Material 

Specific  Heat,  c 

Thermal  Conductivity,  k 

Material  Thickness 

Total  Heat  Transfer  Area 

Frontal-  Area,  A„ 

f  r 
Conductive  Areas  A 
s 
Free  Flow  Area,  A 

Matrix  Flow  Lengtn, 

Matrix  Volume s 

Matrix  Densitys 

Compactness, 

Porosity,  p 

Hydraulic  Diameter 


,A 


m 


m 


m 
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SOLAR  #3 

430  stainless  Steel 

(Btu/lbm  deg  F) 

.11 

(Btu/hr  ft  deg  F) 

12.8 

(in) 

.005 

(sq  ft) 

14.3588 

(sq  ft) 

.06803 

(sq  ft) 

.01325 

(sq  ft) 

.05478 

(ft) 

.24833 

(cu  ft) 

.01689 

(Ibm/cu  ft) 

94.2656 

(aq  ft/cu  ft) 

849.8962 

.80517 

(ft) 

.0037895 

Figure  18.   Geometric  and  Physical  Properties 

For  Solar  #3 
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Matrix  Material 
Specific  Heacs  c 
•   Thermal  Conductivity;,  k 
Material  Thickness 

Total  Heat  Transfer  Area   ,A  (sq  ft) 
Frontal  Area9  Af 
Conductive  Area9rA 
Free  Flew  Area9  A  ? 
Matrix  Flow  Length,  L 
Matr,ix  Volume,  « 
Matrix  Density, 
Compactness a 
Hydraulic  Diameter 


m 
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. SOLAR  #4 

430  stainless  stee 

(Btu/lbm  degF) 

.11 

.(fcfiu/hr  ft 

deg  F) 

12.8 

(in) 

.005 

(sq  ft) 

.14.27694 

(sq  ft) 

.06713 

(sq  ft) 

.01229 

(sq  ft) 

.05484 

(ft) 

.24417 

(cu  ft) 

.01638 

(Ibm/cu  ft) 

88.5918 

(sq  ft/cu  f 

t) 

.81690 

(ft) 

.00375168 

•  Figure  19.   Geometric  and  Physical  Properties 

for  Solar  #4 


60 


m>jK  » 


NaC305R 
5 


SOLAR  #5 

Hatrix  Material 

430  stainless  steel 

Specific  Heat,  c 

(Btu/lbm  deg  F) 

.11 

Thermal  Conductivity,  k 

(Btu/hr  ft  deg  F)  ■ 

12.8 

Material  Thickness 

(in) 

.005 

Total  Heat  Transfer  Area   , 

A    (sq  ft) 

14.4234 

Frontal  Area,  A. 

f  r 
Conductive  Area,  A 

(sq  ft) 

.06781 

(sq  ft) 

'  .01148 

Free  Flow  Area,  A 

(sq  ft) 

.05633 

Matrix  Flow  Lengtfi,  L 

•  (ft) 

.26417 

Matrix  Volume,  V           m 

(eu  ft) 

.0179 

Matrix  Density, 

(lbm/cu  ft) 

81.9345 

Compactness, 

(sq  ft/cu  ft) 

805.13313 

Porosity,  p 

.83Q66 

Hydraulic  Diameter 

.  (ft)   ' 

.0041269 

Figure  20 .   Geometric  and  Physical  Properties 

for  Solar  #5 
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MODIFIED  RECTANGULAR  PASSAGE 


Matrix  Material 
Specific  Heats  e     . 
Thermal  Conductivity 9  k 
Material  Thickness 
Total  Heat  Transfer  Arei 
Frontal  Areas  Af 
Conductive  Area,  A 
Free  Flow  Area9  A 
Matrix  Flow  Leagtn, 
Matrix  Volume,  V  - 
Matrix  Density 
Compactness, 
Porosity,  p 
Hydraulic  Diameter 


L 


m 


(Btu/lbm  deg  F) 
(Btu/hr  ft  deg  F) 
(in) 
,A    (sq  ft) 

(sq  ft)  . 
(s~q  ft) 
(sq  ft) 
(ft) 
(cu  ft) 
(lbm/eu  ft) 
(sq  ft/cu  ft) 

(ft) 


160/40  Perforated  nickel 

.106 
36.0 

.0022 
18.593 

.06953 

.00424 

.0584 

.16098 

.01159 
63.6 
1604.5 

.840 

.001929 


Figure  21.  Geometric  and  Physical  Properties 

for  160/40  T¥ 

Modified  Rectangular 
Passage 
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MODIFIED  RECTANGULAR  PASSAGE 


Matrix  Material 

Specific  Heats  c 

Thermal  Conductivity,  k 

Material  Thickness 

Total  Heat  Transfer  Area, 

Frontal  Area,  A_ 

f  r 
Conductive  Area,  A 

Free  Flow  Area,  A 

Matrix  Flow  Length 

Matrix  Volume,  V 

Matrix  Density  , 

Compactness, 

Porosity,  p 

Hydraulic  Daimeter 


s 


m 


m 


Solid  Nickel 

(Btu/lbm  deg  F) 

.106 

(Btu/hr  ft  deg  F) 

36.0 

(in) 

.002 

(sq  ft) 

20.1875 

(sq  ft) 

.06953 

(sq  ft) 

r     .01014 

(sq  ft) 

.05939 

(ft) 

.16098 

(cu  ft) 

.011588 

(lbm/eu  ft) 

77.7 

(sq  ft/cu  ft) 

1742.1 

.854 

(ft) 

.001961 

Figure  22=,   Geometric  and  Physical  Properties 
for  Solid  Nick«i  Modified 

Rectangular  Passage 
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Figure  23.   Heat  Transfer  and  Flow  Friction 
Characteristics  Solar  #1 
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Figure  24.  Heat  Transfer  and  Flow  Friction 
Characteristics  Solar  #2 
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Figure  25.   Heat  Transfer  and  Flow  Friction 
Characteristics  Solar  #3 
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Figure  26.   Heat  Transfer  and  Flow -Friction 
Characteristics  Solar  #4 
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Figure  27.   Heat  Transfer  and  Flow  Friction 
..  Characteristics  Solar  #5 


68 


tin 


! 

I        !   ;  i:  il! 
Li 


; 


!1!  ! 


MODIFIED     KECT  ..PASSAGE 
160/40  TV  PERFORATED  NICKEL 
.024"  PLATE  SPACING 
POROSITY  .840 

AREA  DENSITY  1604.5 


001 


1000 


5000 


N 


R 


Figure  28.  Heat  Transfer  and  Flow  Friction 
Characteristics  160/40  TV 
Modified  Rectangular 
Passage 
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Figure  29.   Heat  Transfer  and  Flow  Friction 

Characteristics  Solid  Nickel 
Modified  Rectangular  Passage 
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Figure  30.  Comparison  of  Heat  Transfer 

Characteristics  for  all  Surfaces 
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Figure  31.   Comparison  of  Flow  Friction 

Characteristics  for  All  Surfaces 
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Figure  32.   Comparison  of  Flow  Area  Goodness 
Factors  for  All  Surfaces 
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Figure  33.   Heat  Transfer  Power  versus 

Flow  Friction  Power  per  Unit  Area 
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Figure  34.   Heat  Transfer  Data  by"  Means  of 
Cyclic  Technique  Solar  #3 
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Figure  35.   Heat  Transfer  Data  by  Means  of  Single  Blow 
and  Steady  State  Steam  to  Air  Test  Solar  #2 
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HEAT  TRANSFER  and  FLOW  FRICTION  DATA 
SOLAR  #3  BRAZED  TRIANGULAR  FIN 
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Figure  36.   Comparison  of  Experimentally  Determined 
Heat  Transfer  and  Friction  Data  with 
Theory  Solar  #3 
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Figure  37.   Chart  Trace  From  Cyclic  Test 
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Figure  38.      Conduction  Parameter  Analog 
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TABLE  I.      Summary  of  Heat  Transfer  and  Friction  Results 


SCLAR    NC.l    SINUSOIDAL    FINS    SCLIO    NICKEL    5    NILS 


FRICTION    DATA 


HEAT    TRANSFER    DATA 


RUN 


NR 


CONC 
PARAM 


MAX 
SLOPE 


NTU 


NR 


1 

35.66 

.58553 

2 

69.91 

.30134 

* 

1C6.C0 

.20595 

4 

126.80 

.  15680 

5 

170.53 

.  12574 

6 

267.73 

.08029 

7 

394.40 

.05971 

8 

523.69 

.04418 

9 

642.29 

.03663 

10 

912. 84 

.02736 

11 

1164.10 

.02248 

12 

1293.89 

.02064 

13 

1432. 01 

.01873 

14 

1557.13 

.01773 

15 

1678.22 

.01681 

.32719 
.  16642 
.10970 

.08422 
.06764 
.OK312 
.02925 
.022C3 
.01795 
.01263 
.0099C 
.00892 
.008C5 
.00741 
.00668 


.92760 
.93616 
.92300 
.88946 
.86085 
.72688 
.63549 
.56660 
.54626 
.49017 
.39251 
.39769 
.37934 
.37751 
.30071 


25.57 

22.81 

14.74 

11.21 

9.38 

5.48 

2.78 

2.59 

2.09 

1.42 

1.06 

1.C8 

1.03 

1.02 

.00 


35.14 

68.89 

1C4.47 

134.83 

168. oe 

263.87 

388.72 

516.15 

633.05 

899.71 

1147.36 

1275.26 

1411 .41 

1534.71 

1654.05 


.C9482 
.08829 
.05464 
.0*155 
.03475 
.02033 
.01401 
.00959 
.C0775 
.00525 
.00394 
•CC4CC 
.00380 
.00378 
.00000 


FRCNTAL    AREA  AFR 

FLOW    LENGTH  L 

MATRIX    DENSITY  LB/FT2 
FREE    FLOW    AREA  AC 

SCLID    FRCNTAL    AREA  AS 

PCRCSITY  POR 

HEAT    TRANSFER    AREA  A 

HYCRAULIC    RADIUS  RH 

AREA    DENSITY     (A/VM)  (BETA) 

CCNCUCTION    PARAMETER    RATIO    (RS/RK) 


.07022 

.24917 

86.12783 

.05944 

.01078 

.84640 

12.76800 

.00116 

729.70260 

1.0CCC0 


THERMAL    CCNDUCTIVITY  KS 

MATRIX    HEAT    CAPACITY  CS 

ENTRANCE    PRESS. LOSS    CCEF.  KC 

EXIT    PRESSURE       LOSS    COEF.  KE 


36.00 

.1060 

.93767 

-.70431 
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TABLE  II.      Summary  of  Heat  Transfer  and' 'Friction  Results 


SOLAR    NC.2    80/20T    SINUSOIDAL    FINS    SOLID    NICKEL    SPLITTERS    5    MILS 


FRICTION    CATA 


HEAT    TRANSFER    DATA 


RUN 


NR 


CCNC 
PARAM 


MAX 
SLOPE 


NTU 


NR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 


38.67 

75.32 

1  10.01 

139. GC 

212.  84 

276.05 

421.28 

497.17 

570.53 

705.01 

823.83 

975.36 

1096.74 

1368.82 

1538.31 

2317.65 


,53389 
.26494 
.18811 
.  153C4 
.  10594 
.08185 
.05593 
.04840 
.04258 
.03521 
.03072 
.02720 
.02505 
.02142 
.02057 
.01848 


.  14025 
.C72C8 
.04926 
.03968 
.02560 
.01996 
•012S0 
.01096 
.CC963 
.00779 
.00659 
.00562 
.005CO 
.00399 
.00355 
.00226 


.95455 
1  .03187 
.96874 
.91656 
.77382 
.70906 
.58821 
.54974 
.53125 
.50567 
.46394 
.39936 
.38556 
.00000 
.00000 
.00000 


25.60 

21.88 

13.84 

1C.98 

6.60 

5.24 

3.09 

2.35 

1.72 

1.52 

1.30 

1.09 

1.05 

.00 

.00 

.00 


38.12 

74.24 

109.21 

136.98 

209.78 

272.04 

415.24 

490.02 

562.26 

694.80 

812.02 

961.25 

1080.68 

1349.10 

1516.14 

2284.25 


.10087 
.08623 
.05455 
.04330 
.02602 
.02066 
.01217 
.00926 
.00679 
.00598 
.00513 
.00428 
.00412 
.00000 
.00000 
•00000 


FRCNTAL    AREA 

FLCW    LENGTH 

MATRIX  DENSITY 

FREE  FLOW  AREA 

SCLID  FRCNTAL  AREA 

PCROSITY 

HEAT  TRANSFER  AREA 

HYDRAULIC  RADIUS 

AREA  DENSITY  (A/VM) 

CONDUCTION  PARAMETER 


AFR 
L 
LB/FT3 
AC 
AS 
POR 
A 

RH 
(BET/) 
RATIO  (RS/RK) 


.07C22 
.25583 

70.30855 
.05989 
.01C34 
.85281 

12.093C2 

.00116 

673.13967 

.46730 


THERMAL  CCNCUCTIVITY  KS 

MATRIX  HEAT  CAPACITY  CS 

ENTRANCE  PRESS. LCSS  CCEF.  KC 

EXIT  PRESSLRE   LCSS  CCEF.  KE 


36.00 

•  1060 

.07222 

.02166 
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.TABLE  III.      Summary  of  Heat  Transfer  and  Friction  Results 


SCLAR    NO. 2    SINUSOIDAL    FINS    SOLID    430    STAINLESS    NAX    eRAZE    5    NILS 


FRICTION    DATA 


HEAT    TRANSFER    DATA 


RUN 


NR 


CONC 
PARAM 


MAX 
SLOPE 


1 

60, 

.19 

.20320 

2 

<n , 

.18 

.  13761 

3 

121. 

.78 

.  10327 

4 

172. 

.54 

.07288 

5 

221  , 

.36 

.06223 

6 

342, 

.40 

.04132 

7 

456, 

.46 

.03125 

8 

567, 

.42 

.02574 

9 

672. 

.09 

.02200 

10 

902, 

.44 

.01737 

11 

1194 

.66 

.01378 

12 

1436, 

.51 

.01193 

13 

1921, 

.60 

.00975 

FRONTAL    AREA 
FLOW    LENGTH 
NATRIX    DENSITY 
FREE    FLOfc    AREA 
SOLID    FRONTAL    AREA 
POROSITY 

HEAT    TRANSFER    AREA 
HYDRAULIC    RADIUS 
AREA    DENSITY    (A/VM) 


.07428    1.03852 
.049C7       .97702 


.03673 
.02593 
.02021 
.013C8 
.0098C 
.00789 
.006^6 
.00496 
.00375 
.00312 
.00233 


AFR 
L 
LB/FT2 
AC 
AS 
PCR 
A 

RH 
(BETA) 


.89137 
.77241 
.68751 
.57362 
.52398 
.49362 
.45011 
.00000 
.00000 
•CCOOO 
.00000 


NTU 


23.20 
14.24 
9.98 
6.57 
4.82 
2.81 
1.65 
1.44 
1.25 
.CO 

•  CO 
.00 

•  CO 


CCNOUCTICN  PARAMETER  RATIO  (RS/RK) 

THERMAL  CONDUCTIVITY  •  KS 

NATRIX  HEAT  CAPACITY  CS 

ENTRANCE  PRESS. LOSS  CCEF.  KC 

EXIT  PRESSURE   LOSS  CCEF.  KE 


.06803 

.24833 

94.26563 

.05478 

.01325 

.8C517 

1U.35884 

.0CC95 

849.89621 

1.0G000 

12.80 

.1100 

.97439 

-.65449 


NR 


59.22 

89.87 

120.03 

170. C6 

218.18 

337.48 

449.90 

559.26 

662.44 

889.47 

1  177.47 

1415.83 

1893.93 


.07076 
.04325 
.03021 
.01995 
.01464 
.00853 
.00502 
.00438 
.00379 
.CCCCO 
.OOCCO 
.COOCO 
.OOCCO 
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TABLE  IV.      Summary  of  Heat  Transfer  and  Friction  Results 


SOLAR    NC.M    SINUSOIDAL    FINS    SOLID    430    STAINLESS    PATL.    5    MILS 


FRICTION    DATA 


HEAT    TRANSFER    DATA 


RUN 


NP 


CONC 
PARAM 


MAX 
SLOPE 


NTU 


NR 


1 

30.98 

.37693 

2 

59.13 

.20381 

3 

88.76 

.  14272 

u 

1  12.87 

.  11402 

5 

171.76 

.07412 

6 

215.99 

.05929 

7 

338.77 

.03983 

8 

447.77 

.03100 

9 

557.45 

.02546 

1C 

664.38 

.02172 

11 

888.19 

.01741 

12 

1162.42 

.01426 

13 

1415. 3M 

.01239 

14 

1908.14 

.01161 

.  1341)5 
.CICUk 
-.0U6S2 
.U36E9 
.02425 
.01929 
.01220 
.00929 
.007H7 
.C0627 
.00469 
.00358 
.00296 
.00219 


1  .16696 
1.07385 
.97752 
.89210 
.74849 
.68279 
.5641  1 
.5078C 
.48585 
.41536 
.00000 
.00000 
.00000 
.00000 


2.00 

27.52 

14.10 

1C.01 

6.04 

4.74 

2.64 

1.53 

1.4G 

1.13 

.00 

.00 

.00 

•  CO 


30.53 

58.28 

87.49 

111.25 

169.30 

212.89 

333.91 

441.35 

549.45 

654.85 

875.45 

1145.75 

1394.95 

1880.66 


.00611 
.08412 
.04310 
.03061 
.01645 
.01449 
.CC8C6 
.00467 
.00429 
.00346 
•COCCO 
.000C0 
.OOOCC 
.00000 


FRCNTAL    APEA  AFR 

FLOW    LENGTH  L 

MATRIX    DENSITY  LB/FT2 

FREE    FLOW    AREA  AC 

SCLID    FRCNTAL    AREA  AS 

PCROSITY  PCR 

HEAT    TRANSFER    AREA  A 

HYDRAULIC    RADIUS  RH 

AREA    DENSITY    (A/VM)  (BETA) 
CCNDUCTICN    PARAMETER    RATIO    (RS/RK) 


.06713 

.24417 

88.59182 

.05484 

.01229 

.81690 

14.27694 

.00094 

870.99223 

l.OCOOO 


THERMAL    CCNCUCTIVITY  KS 

MATRIX    HEAT    CAPACITY  CS 

ENTRANCE    PRESS. LCSS    CCEF.    KC 
EXIT    PRESSURE       LCSS    CCEF.    KE 


12.80 

.1100 

.96374 

-.66901 
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TABLE  V. 


Summary  of  Heat  Transfer  and  Friction  Results 


SCLAR    NC.5    TRIANGULAR    FINS    SOLID    430    STAINLESS    NATL.    5    MILS 


FRICTION    DATA 


HEAT    TRANSFER    OATA 


RUN 


NR 


CONC 
PARAM 


MAX 
SLOPE 


NTU 


NR 


1 

63.05 

.22228 

2 

96.54 

.  14597 

3 

121.54 

.  11431 

y 

183. 46 

.08385 

5 

231.1  1 

.06594 

6 

360.72 

.04488 

7 

482.55 

.03439 

8 

600.21 

.02846 

9 

708.74 

.02472 

1Q 

964.34 

.01907 

11 

1247.65 

.01548 

12 

1548.85 

.01341 

13 

2040.53 

.01207 

.061  ie 

.03995 
.03172 
.021C1 
.01667 
.01068 
.00798 
.00642 
.00544 
.004C1 
.C031C 
.00250 
.00190 


1.04806 
.95656 
.88205 
.73707 
.66779 
.56045 
.50007 
.47897 
.45696 
.00000 
.00000 
.00000 
.00000 


21.15 

12.62 

9.54 

5.81 

4.4e 

2.57 

1.48 

1.37 

1.28 

.00 

.00 

.00 

.00 


62.15 

95.15 

119.79 

180. e3 

227.80 

355.54 

475.63 

591.59 

698.57 

950.46 

1229.68 

1526.55 

2011.13 


.06573 
.03924 
.02967 
.01805 
.01393 
.008C0 
.00460 
.00426 
.00396 

.oooco 
.oooco 

.00000 
•00000 


FRONTAL  AREA  AFR 

FLCW  LENGTH  L 

NATRIX  DENSITY  .  LB/FT2 

FREE  FLOfo  AREA  AC 

SCLID  FRCNTAL  AREA  AS 

PCROSITY  PCR 

HEAT  TRANSFER  AREA  A 

HYDRAULIC  RADIUS  RH 

AREA  DENSITY  (A/VM)  (BET/) 

CCNDUCTICN  PARAMETER  PATIO  (RS/RK) 


.06781 

.26417 

81.93452 

.05633 

.01148 

.83066 

14.42340 

.0C1C3 

805.13313 

1.CCCC0 


THERMAL  CCNDUCTIVITY  KS 

MATRIX  HEAT  CAPACITY  CS 

ENTRANCE  PRESS. LOSS  CCEF.  KC 

EXIT  PRESSURE   LOSS  CCEF.  KE 


12.80 

.1100 

.95140 

-.68569 
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TABLE  VI*      Summary  of  Heat  Transfer  and  Friction  Results 


160/40TV    PERFORAiED    NICKEL    FINS    AND    SPLITTERS 


FRICTION  DA  i A 


HEAT  TRANSFER  DA1A 


RUN 


NR 


CCND 
PARAM 


MAX 

SLOPE 


NlU 


NR 


812.99  .02545 


2 
3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 


653 

527 

406 

310 

206 

124 

69 

39 

69 

58 

48 


01 
22 
11 
00 
38 
11 
91 
33 
29 
63 
89 


03104 
03751 
04445 
06028 
08782 
13378 
23837 
42495 
23339 
27522 
34594 


38.50  .41329 


70031a"  .55590    2.54    801.26  .00635 

00816 
01091 
01287 
01658 
02405 
03932 
06619 
00500 
07539 
00500 
00500 
37.94  .00500 


00395 

.59357 

3. 

27 

643 

.60 

00490 

.65507 

4. 

,36 

519 

.61 

00635 

.69880 

5. 

,  15  " 

400 ; 

.27 

00832 

.77426 

6. 

,63 

305, 

.54 

01253 

1 

.89948 

9. 

62 
,73 

203, 
122 

.40 

02080 

.06826 

15, 

.32 

03694 

1 

.17903 

26, 

,47 

68 

.91 

06561 

1 

.16650 

2, 

,00 

38, 

.76 

03736 

1 

.21282 

30, 

15 

68 

.29 

04414 

1 

.23549 

2, 

,00 

57 

.78 

05293 

1 

.22826 

2. 

,00 

48, 

.19 

06718  1.19438    2.00 


FRONTAL  AREA 
FLOW  LENGTH 


MATRIX  DENSITY 

FREE  FLOW  AREA 

SOLID  FRONTAL  AREA 

POROSITY 

HEAT  TRANSFER  AREA 

HYDRAULIC  RADIUS 


AFR 

L 


.06264 
.16670 


(  L  8  /  F  T  3  ) 
AC 
AS 


RCR 
"A" 

Rh 


AREA  DENSITY  (A/VM) 

COND. PARAMETER  CORRECTION 


SETA 
t 


70.51272 

"TO  58  4  0" 

.0  0424 

.84000 
18.59300 

.00048 


1604  *  5 
.82844 
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TABLE  VII.      Summary  of  Heat  Transfer     and  Friction  Results 


002    IN.     SOLID    NICKEL    MOD. PARALLEL    PLATE 


FRICTION    DAiA  HEAT    TRANSFER    DA1A 


CGND      MAX         

RUN      NR      F  PARAM    SLOPE    RTTJT      NR 


1  740.42  .02836        .00925  .57103  2.79    729.78  .00653 

2  608.95  .03389        .01125  .58706  3.08    600.20  .00722 

3  481.75  .04027 .01420  .61318    3-52 474.84  .00825 

4  373.92  .04731"       .01832  .67679  4.63"   368.55  .01085 

5  286-81  .05896        .02388  .71948  5.43    282.69  .01271 

6  187.14  .08669 .03660  .81839    7.74 184.45  .01812 

7  115.43  .12779"       .05931  . 86112  9 . 27    113.77  .02170 

8  65.08  .21276        .10518  .81118  7.78     64.15  .01823 

9  35.72  .37233 119163  .77080    5.55 35.21  .01299 

FRONTAL  AREA               AFR  .06953 


FLOW  LENGTH  L  .16670 

MATRIX  DENSITY  fLtJ/Fl3)       77763067' 

FREE  FLOW  AREA  AC  .05939 


SOLID  FRONTAL  AREA         AS  .01014 

POROSITY        ^CR  .85420 

HEAT  TRANSFER  AREA         A  20.18750 


HYDRAULIC  RADIUS  RH  .00049 

AREA  DENSITr  (A/VM)        BETA  1742.1 

COND. PARAMETER  CORRECTION  1.00000 
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APPENDIX  A 

DATA  REDUCTION  RELATIONSHIPS 

GEOMETRY 

Porosity  ■  Matrix  free  flow  area 
Matrix  frontal  area 


Af> 

Hydraulic  diameter  ■  4x  Free  flow  area 

Wetted  perimeter 

4Ac_L 

Area  density  ■  Heat  transfer  surface  area 

Matrix  volume 


Ac 


DH  -  *JTT*Z  (A-2) 


G  A       A 

MASS  RATE  OF  FLUID  FLOW 

The  mass  flow  rate,  m,  is  calculated  in  accordance  with  references 
[21]  and  [23]. 

m     -  359  k  doZ  Fc,    Y  VAP,  V     Obw/hf)  (A-4) 

Where: 

K  =  J.  _jol      the  flow  coefficient  including  velocity  of  approach 

C  =  orifice  coefficient  of  discharge  from  [21] 

]P  =  ratio  of  orifice  diameter  to  inside  pipe  diameter 

C«o=  orifice  diameter 

U,,=  thermal  expansion  factor 

Y  ■  expansion  factor 

p 
\    -       «tp    the  specific  weight  of  the  fluid  flowing 

P  =  absolute  static  pressure  at  the  orifice  (lbf/sq  ft) 

R  =  53.35  ft  lbF,  the  universal  gas  constant  for  air 

lbm°R 
T  =  absolute  temperature  (°R)  at  the  orifice 
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Ar0   m       Pressure  drop  across  the  orifice  in  inches  of  water 
Making  substitutions  for  s.     and   [   ,  (A-4)  becomes: 

=     4F?& VAP°  RT  (A-5) 

In  accordance  with    [23]    figures    (38)   and   (40b): 
Fft.ftJl.0 

Also: 

*  -  (  Patti  -  rf#6   )  -^12  X  144       (lbf/sq  ft) 
Where:        R^th  is  measured  in  inches  of  mercury 
Po   is  measured  in  inches  of  water 

and: 

T  -  to     +  459.7 

where  ^0  is  measured  in  degrees  Fahrenheit 

Introducing  the  values  and  expressions  above  into  (A-5)  produces 
the  working  formula: 

^  s  ?*V'  **%/**  fe™l&)        (A-e: 


REYNOLD'S  NUMBER 


By  definition: 


Subs 


tituting      4AcL      for  DH    »   and    ^-    for    ©        I 


(A-7) 
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FANNING  FRICTION  FACTOR 


EZZZZZZZZZZZZZZZZ22I 

mnmznzznzzznm 


mznzzzzznzzzzzzm  . — rv 

a   VIM  gg  Zg  "/WW  *Wi  b  Z      I 1/ 


The   following  equation  describes   the   flow  through  th#  matrix    \\1\\ 

_     EMT«AWCff  ACCEi-fRATlOfJ  F*.ICTfO»J  £xiT 

where  p  =  porosity  (   /^f  ) 

It  is  assumed  that  'U&.^V,  and  'V^'x.  IT^  since  the  associated 

pressure  changes  are  very  small  compared  with  the  total  pressure.   It 

DT 
is  also  assumed  that  the  perfect  gas  low  applies  (  1T=  -=-  ).  Making 

the  necessary  substitutions  and  solving  for   f  produces  for  the  iso- 
thermal case: 

f -  [«*.  #  -  wx^teHWMH1?  <a-8) 


Since  the  first  term  within  the  brackets  above  is  much  larger  than 
the  other  terms,  the  following  approximation  can  be  made: 
Pi* Pi  _  p   ^  p  ~  "P 

Thus  equation  (A-8)  becomes: 


f« 


**•««  §  -(«c  +  /<e)  +  ^  (  1  +  P1)]  *f     (A-9) 
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k  c   and  /<tt  ,  the  entrance  and  exit  loss  coefficients  are  functions  of 
porosity,  matrix  configuration,  and  Reynold's  number,  and  are  obtained 
in  accordance  with  reference  [11]. 

MAXIMUM  SLOPE 

The  maximum  slope  of  the  generalized  cooling  curve  as  previously 
mentioned,  is  a  function  of  Ntu  and  conduction  parameter,  X  : 

MX 


(A-10) 


where: 


>M-5 


WsC 


s^-S 


Therefore: 


and    (s/tH  =  M. 


and: 


fA-n) 


also: 


and: 


***-£*    .   tlFu-  fcpj 


tf,  -t; 


combining  (A-ll)  and  (A-12)  yields: 

V/*JL*2 


+  1 


-l^<M.-*rfj 


-  ^s 


dft*- 


Cf  tp.-fc; 


max 


'A 


(A-ia) 


(A-13) 
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Now  consider  the  following  sketch  of  a  generalized  cooling  curve: 


t 


e 


d(tFl-tFl) 


d© 


HAX 


y. 


69 

d(tPt-tFl) 
(tr,"ti) 


CHAftT  3P£E& 

=   Y 

=  m 


Combining  the  above  produces  the  maximum  slope: 


=  —  i-t-i-Z-tCMILT  SPEED 
MAK  (A-14) 


Ntu  can  now  be  determined  from  conduction  parameter  X  and  maximum 
slope. 
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COLBURN   j 

Colburn's  heat  transfer  parameter,  j,  is  defined: 

J  ■  NM  NPrVa 


'*■  (A-15) 


thus: 


j-Nfcu^Nr/'3  <A-i6) 


From  equations  (A-l)  and  (A-3):      ^c        J3         A 

j«Mt«MpfV»^-|  (A-17) 


From  equation  (A-17)  it  can  be  seen: 

HEAT  TRANSFER  POWER 

The  heat  transfer  power  per  unit  area  per  degree  temperature 
difference  is  [12]: 

h  =    TZ*  Oh   J  Nr  <a"18) 


For  dry  air  at  standard  conditions  equation  (A-18)  becomes: 


U«.«rihi.      Bt^<trt^r)    (A,9) 
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FLOW  FRICTION  POWER 


The  flow  friction  power  per  unit  area  is  defined  as  [12]: 
For  dry  air  at  standard  conditions  equation  (A-20)  becomes: 

ESTo='-»"<0-7fo)3f(£§o)3      %„      u-n, 


For  comparison  purposes  the  surface  geometries  were  reduced  to  a 
common  hydraulic  diameter. 
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APPENDIX  B 
DESCRIPTION  OF  EQUIPMENT 


PRIME  MOVER: 


Air  as  the  working  fluid  is  drawn  through  the  system  by  means  of 
a  30  HP,  multistage,  Spencer  Turbo-Compressor.   This  compressor  is  de- 
signed for  a  flow  rate  of  550  cubic  feet  per  minute. 

FLOW  METER: 

Flow  rates  are  measured  by  means  of  an  ASME  standard  orifice  section 
designed  for  convenient  changing  of  orifice  plates.   Thin  plate  concen- 
tric orifices  with  throat  diameters  ranging  from  .308  to  2.31  inches  are 
employed  allowing  for  a  wide  range  of  flow  rates.   Pressure  taps  are 
located  d  upstream  and  d/2  downstream.   Inside  pipe  diameter,  d,  is 
3.08  inches. 

HEATER  SYSTEM: 

The  heater  section  is  made  up  of  28  grid  type  nichrome  wire  heaters. 
Two  heaters  are  wound  to  each  bakelite  frame  and  thus  may  be  energized 
in  pairs.  A  schematic  wiring  diagram  for  one  heater  frame  is  shown  below: 


^ 


HEA7E*.  SVSTcM 
0*i-   OFF 

SWiTC/f 


VOLTAGE 
SOURCE 


L 


HEATetZ 
SWtTfcH 


1. 

I*/ 

1 


WEXT      -. 
U£*TZQ3 


S 
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When  energized,  all  heaters  operate  in  parallel,  the  number  of 
heaters  in  use  depending  upon  flow  rate. 

For  maximum  slope  operation  the  voltage  source  is  220  volts  AC  and 
can  be  varied  by  means  of  a  General  Radio  Company  "Variac"  Autotrans- 
former. 

For  cyclic  operation  the  voltage  source  is  a  1.5  KW  General  Electric 
Amplidyne  connected  in  series  with  a  D.C.  source.   The  resulting  input 
to  the  heaters  is  a  sine  wave  voltage  with  an  amplitude  of  50  volts  and 
an  average  value  of  50  volts.   Frequencies  of  .1  to  1.0  cycles  per  second 
are  employed  with  a  much  wider  range  available. 

MATRIX  HOLDER  AND  TEST  SECTION; 

This  section  is  constructed  of  polyethylene  machined  to  close  toler- 
ances.  The  T_  and  T,  thermocouple  grids  and  the  matrix  upstream  and 
downstream  pressure  taps  are  mounted  in  the  frame.   The  T-  thermocouple 
grid  is  mounted  in  the  removable  matrix  holder.   %"  sheet  styrofoam  in- 
serts provide  insulation  and  a  snug  fit  for  the  test  matrix.   The  flow 
channel  measures  3-3/16"  X  3-3/16"  and  matrices  up  to  3"  in  length  may 
be  accommodated  in  the  test  section.  (See  Figure  12) 

TEMPERATURE  MEASUREMENT: 

The  temperature  at  the  orifice  meter  is  measured  by  means  of  a 
copper-constantan  thermocouple  and  read  in  millivolts  with  a  Rubicon 
Portable  Precision  Potentiometer. 

All  other  temperatures  are  measured  by  means  of  iron-constantan 
thermocouples,  each  temperature  being  measured  by  a  group  of  five 
thermocouples  connected  in  series.   The  T.  thermocouples  are  bound  to- 
gether and  insulated  from  each  other  by  teflon  tape.   This  bundle  is 
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mounted  at  the  throat  of  the  intake  cone,  upstream  of  the  wire  screen 
flow  straightener.   The  T_,  T,-,and  T  grids  are  connected  to  a  multiple 
barrel  switch  so  that  the  emf  generated  by  T1  may  oppose  the  emf  generat- 
ed by  either  T.  or  T~.   Thus  the  difference  between  ambient  and  matrix 
inlet  or  exit  temperature  can  be  read.   These  temperature  differentials 
are  recorded  on  a  Minneapolis-Honeywell  MBrown"  Strip  Chart  Recorder  with 
variable  sensitivity  and  chart  speed.   For  maximum  slope  operation,  T. 
versus  T_  is  recorded.   For  cyclic  temperature  variation,  T.  versus  T~, 
and  T.  versus  T~  are  recorded. 

PRESSURE  MEASUREMENT; 

The  four  pressure  taps  shown  in  Figure  8  are  connected  to  Ellison 
draft  gages  by  means  of  Imperial  Company  "Poly-flo"  plastic  tubing.   The 
draft  gages  are  a  combination  of  inclined  and  vertical  gages.   They  are 
shown  in  Figure  9.   For  small  pressure  differentials,  Vernon  Hill  Company 
Type  "C"  Micromanometers  are  used.   Pressures  of  1.25"  to  3.0"  of  water 
are  read  by  means  of  the  inclined  draft  gages.  Pressures  below  this  range 
are  measured  by  means  of  the  micromanometers;  Pressure  above  3.0"  are 
read  on  vertical  draft  gages. 
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APPENDIX  C 
CONDUCTION  PARAMETER  FOR  PERFORATED  MATERIAL 
The  porous  material  used  in  the  perforated  plate  matrices  investi- 
gated in  this  report  was  made  by  Perforated  Products,  Inc.  by  means  of 
an  electro-depositing  process. 

For  a  solid  plate  or  fin,  the  conduction  parameter  is  defined  as: 


_  fcs/4s 


LC 

where  A  is  the  solid  cross-sectional  area  for  conduction  which  is 
s 

identical  to  the  product  of  plate  width  and  thickness,  and  L  is  the 
conduction  path  length  which  is  identical  to  the  matrix  flow  length. 

For  a  perforated  material  such  as  the  type  160/40  TV  shown  in 
Figure  38(a),  the  cross-sectional  area  for  conduction  varies  in  the 
flow  direction,  and  the  conduction  path  length  is  greater  than  the  flow 
length. 

The  definition  of  conduction  parameter  is  based  upon  a  solid  materi- 

a1'  i*e"        ^    _i_ 
*~  LC     '    KykC 

where  R^,  ■  -. — r    the  matrix  thermal  conduction  resistance, 
th   k  A 
s  s 

Using  the  above  definition  for  a  perforated  material,  an  "equivalent" 
conduction  parameter,  Xfc  ,  can  be  evaluated  by  determining  the  thermal 
resistance  (R^v),  .  By  defining  a  relative  shape  factor  S1,  the  conduc- 
tion parameter  of  a  perforated  material,  Xfc  ,  can  be  related  to  the  con- 
duction parameter  of  the  same  material  as  if  it  were  unperforated: 
Xk  ■  ^'X  *  *• 


a  j<  ■  >-»  /\    *• 
or: 
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By  the  analogy  between  conduction  heat  flow  and  electric  current,  note 
that  for  similar  shapes,  the  ratio  of  electrical  resistances: 

£  -    fijfa.     «  s' 

In  Figure  38(b)  are  shown  two  shapes,  each  representing  the  shaded 
portion  of  Figure  38(a).   The  shape  on  the  left  represents  an  unperfor- 
ated  material,  and  the  shape  on  the  right  is  an  incremental  model  for  the 
160/40  TV  material.   The  relative  shape  factor  is  determined  from  the 
ratio  of  the  electrical  resistances  of  these  two  models.   Since  the 
material  is  of  uniform  thickness,  two  dimensional  models  can  be  construct- 
ed from  electrically  conducting  paper  (such  as  Teledeltos  paper)  as  shown 
in  Figure  3S*tc) ,  with  low  resistance  silver  paint  applied  as  indicated. 
The  electrical  resistances  are  then  measured  with  an  ohm-meter.   This 
analog  method  produced  the  following  results  for  the  type  160/40  TV  per- 
forated nickel  material: 

-'•JBu  *  £  «  ££=  .3/5" 

thus:     Xk=  .ZIS"\ 

Whereas,   Nk  =  .127^    as  determined  in  reference  [22]  using: 

where:     .     _.  . 

and:    A^  and  k  .  are  approximations  indicated  in  Figures  38(a),  and 
(b). 

The  effects  of  these  approximations  upon  Ntu  and  j  are  shown  in 
Figure  15. 
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APPENDIX  D 
EXPERIMENTAL  UNCERTAINTY  INTERVALS 

According  to  reference  [13],  the  uncertainty  interval  for  a 
result  is: 

*-[«;^ffliwo** *fifc*y]*   <-> 

w-i »  w  , . .  .w  are  the  uncertainty  intervals  of  the  factors  making 
up  the  result. 

R  is  the  expression  for  the  result 

V.,  V_,...V  are  the  factors  making  up  the  result. 

l   z     n 

The  expression  for  Colburn  j  is: 

The  uncertainty  interval  for  Colburn  j  is: 

y 

*-  [  (Ifc  w-)%  (St^T  +<8W  *gfcft  w^)z] * 


For  Ntu  -  3.0;    A=  0  etc: 
^i     =  7.57, 

For  Ntu   -   25.0;    X«   .05  etc: 

^   »    10.2% 
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APPENDIX  E 
DIGITAL  COMPUTER  PROGRAM  FOR  DATA  REDUCTION 

Included  herein  is  a  computer  program  for  a  CDC  1604  digital  com- 
puter.  By  means  of  this  program  all  necessary  calculations  were  made 
to  produce  the  heat  transfer  and  flow  friction  data  appearing  in  the 
tables  and  graphs  of  this  report. 

A  data  reduction  program  was  written  by  Bannon  [1],  This  program 
was  then  further  sophisticated  by  Piersall  [22],  producing  conduction 
parameter  \  ,  maximum  slope,  and  friction  factor  as  well  as  other  flow 
information  from  geometric  and  experimental  data  inputs. 

The  present  program  will  now  produce  all  necessary  heat  transfer  and 
flow  friction  data  in  table  form  from  a  minimum  of  input.  When  data  from 
the  single-blow  technique  is  processed,  a  curve- fitting  interpolation 
subroutine  is  utilized  to  determine  Ntu  from  maximum  slope  and  conduction 
parameter.  The  input  for  this  subroutine  is  Howard's  f9]  maximum  slope 
data  (Table  VIII).  For  the  cyclic  technique,  Ntu  is  determined  from 
equation  (2-33)  by  means  of  an  iteration  process. 

The  program  will  process  data  obtained  by  means  of  either  the 
single-blow  or  cyclic  technique  and  data  from  more  than  one  matrix  can  be 
reduced  in  each  computer  run. 
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figuration.  The  modified  rectangular  passage  surfaces  showed  better  overall 
performance  than  the  triangular  fin  surfaces,  and  matrices  constructed  from 
perforated  material  showed  better  overall  performance  than  similar  surfaces 
constructed  from  non-perforated  material.   For  similar  matric  configurations, 
small  variations  in  manufacturing  procedure  had  little  effect  upon  heat 
transfer  and  flow  friction  performance. 

The  heat  transfer  data  was  obtained  by  means  of  the  single-blow  transient 
testing  technique.   A  cyl.ic  method  of  transient  testing  was  also  investigated, 
and  found  to  be  both  reliable  and  useful  in  that  the  additional  method  of 
testing  extends  the  Reynolds  number  range  for  the  experimental  apparatus. 
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